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We study precursors to the global market crash that occurred on all main stock exchanges 
throughout the world in October 2008 about three weeks after the bankruptcy of Lehman 
Brothers Holdings Inc. on 15 September. We examine the collective behavior of stock returns 
and analyze the market mode, which is a market-wide collective mode, with constituent 
issues of the FTSE 100 index listed on the London Stock Exchange. Before the market 
crash, a sharp rise in a measure of the collective behavior was observed. It was shown to 
be associated with news including the words "financial crisis." They did not impact stock 
prices severely alone, but they exacerbated the pessimistic mood that prevailed among stock 
market participants. Such news increased after the Lehman shock preceding the market 
crash. The variance increased along with the cumulative amount of news according to a 
power law. 



§1. Introduction 



What causes a drastic price change such as that occurring during a crash? Ac- 
cording to the efficient market hypothesis, a prevailing paradigm of mainstream 
economics, stock prices change because of news that comes as a surprise to mar- 
ket participants. However, several previous works suggest that this picture of price 
changes is not empirically acceptableSil 

From the Paribas' shock in August 2007 until the beginning of 2009, many stock 
markets throughout the world experienced considerable stock price declines because 
of ripple effects of the US sub-prime loan problem and the subsequent financial crisis. 
In Figure Qa) , the time series of the FTSE100 index and the daily log-return are 
shown for the period from May 2007 through January 2009. Regarding some heavy 
falls such as that of the Paribas' shock on 9 Aug 2007, the buyout of Bear Stearns 
by JP Morgan on 17 March 2008, and the Lehman shock on 15 September 2008, 
they were probably caused by large exogenous shocks. However, larger falls such as 
the crashes of January and October 2008 are related to no specific news that might 
equal the magnitude of the subsequent drop-off. As an example, the evolution of 
the 1-min log-return of Royal Dutch Shell pic (RDSA) and the daily frequency of 
its large values of amplitudes larger than the twofold standard deviation during the 
period are also shown in Figure [l|bf^| It is apparent that the intermittent large 
price changes occurring during the daytime around 21 Jan. and 8 Oct. 2008 market 
crash without specific news are more frequent and persistent compared to the days 
showing steep declines that are apparently caused by important economic news. 

Especially for the market crash of Oct. 2008, the daily frequency of large returns 
gradually grows from the Lehman shock on 15 Sep. until the maximum on 8 Oct. 
and decays slowly. Results demonstrate that the relaxation dynamics of a financial 



*' Here, we exclude overnight price changes. 



typeset using -pTpf^ ds (Ver.0.9) 



2 



J. Maskawa 





Fig. 1. Large price declines during May 2007 to Jan 2009. (a) Time series of the FTSE100 index 
(black line) and the daily log-return (gray line) . (b) Evolution of the 1-min. log-return of Royal 
Dutch Shell pic (RDSA) (total: 198490 min)(gray lines) and the daily frequency of large values 
with amplitudes greater than the twofold standard deviation during the period (black lines). 
Overnight price changes are excluded. 



market immediately after the occurrence of a crash resemble an earthquake after- 
shock® The frequency of a large aftershock decays according to a power law. The 
power law of earthquake aftershocks was reported by Omori in 1984Hf The market 
behavior before and after the scheduled macro-economic news announcements, such 
as U.S. Federal interest rate changes, has also been studied, revealing a similar law 
of foreshoc k© Figure k| portrays a plot of the cumulative frequency N(t) of large 
log-returns with amplitudes larger than a given threshold against the elapsed time 
t, which has an origin coordinate that is the time of the largest amplitude of the 
log-return during the daytime of 8 Oct. 

According to previous work by Petersen et al.jlf the cumulative frequency N(t) 
is fitted by the following power function: 

\N(t) -N(0)\ ~ \tf> (i = b,a), (1-1) 
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Fig. 2. Cumulative frequency of large 1-min. log-returns of RDSA. From top to bottom, we respec- 
tively show cumulative frequencies for the threshold 4ctm, 5a m, 6&m, 7&m (black solid lines). 
Here om is the standard deviation during May 2007 - Jan 2009. The dashed lines are the best 
power-law fits by Eq. ( | 1 ■ 1[ | . The numbers in parentheses attached to each curve are {[3b, Pa) for 
the corresponding threshold. The time-series of 1-min. log-return is also shown (gray line). 



where the time regions t < (i = b) and t > (i = a) respectively correspond to 
a foreshock and aftershock. We have the inequality of < fib < fi a < 1, which is 
the same result as that for the scheduled macro-economic newso Non-vanishing fib 
means that a foreshock occurs before the main shock. 

Can we predict a singular point where the curvature of the power function change 
discontinuously at the main shock? In principle, it is impossible to determine the 
singular point because the estimated value fib would not be the value if the time of 
the main shock were not there. In this paper, we examine the collective behavior 
that affects the stock prices. We analyze the time-series of the average of normalized 
stock returns of stock index portfolio, which is approximately the first principal 
component of multivariate time-series of those returns and represents the market- 
wide collective behavior of stock prices0 Designated as the "market mode" here, 
its rigorous definition is provided in the next section. 

As described in this paper, using the framework of the Multifractal Random 
Walk (MRW) model|5 we show the upwelling of collective behavior before the market 
crash. That is related to the news which includes the word "financial crisis", which 
has no severe impact to stock price alone but which contributes to worsening of 
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the pessimistic mood among stock market participants. The amount of such news 
increased after the Lehman shock, which preceded the market crash. The variance 
increases along the cumulative amount of news according to a power law. 



§2. Data 



We analyze the multivariate time-series of stock returns of 105 selected issues 
from the constituent of FTSE100 index listed on the London Stock Exchange for the 
period from May 2007 to January 2009 (445 days j^J which includes the period during 
which drastic price changes occurred because of the US sub-prime loan problem and 
the subsequent financial crisis. In Figure [T^a), the time series of the FTSE100 index 
and the daily log-returns are shown for the period. We exclude the ove rnight price 



changes and specifically examine the intraday evolutions of returns**) The total 
length of the time-series investigated here is T = 198, 490 min. The return of the 
issue i is the difference of the logarithms of log-price of t and t — 5t: 

SXi(t) = log(Pi(t)) - log(Pi(t - St)). (2-1) 

Here, we set 5t = 1 and use the 1-min log-return. Next, we define the market mode 
M{t) as the average of normalized returns: 

^) = ^E— > (2-2) 

where the number of issues N = 105 and is the standard deviation of log-return 
of i-th issue during the period. The market mode 6M(t) is approximately the first 
principal component of multivariate time-series of returns and represents the market- 
wide movement of stock prices Finally, to remove the effect of the intraday U shape 
pattern of market activity from the time-series of the market mode, 5M(t) is divided 
by the standard deviation of the corresponding time of day, which is obtained from 
the whole of the time-series. We use the same notation M(t) for the de-seasonalized 
market mode. 



§3. Empirical results 

We use the framework of the Multifractal Random Walk (MRW) mo de§ to an- 
alyze the time-series of the de-seasonalized market mode M(t). We briefly introduce 
the MRW model in Appendix A, and show that the market mode is a multifractal 
random process and that it is well fitted by the model. The market mode is ap- 
proximately the first principal component of multivariate time-series of N returns 
and represents the market-wide collective behavior of stock prices In the principal 

*' The constituents of FTSE100 index are updated frequently. We select those issues which had 
appeared on the list of the constituents at least once during the period and which had been listed 
on the London Stock Exchange throughout the period. 

**' To average out the normalized returns of selected issues, we skip 60-110 min immediately 
after the market opening, varying from day to day until all the issues get opening prices. 
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component analysis, the variance of principal components expresses the contribution 
to all the variation of multivariate time-series. Prom equations (A-5) and (A-6), the 
variance of ugt the logarithm of the volatility is expressed by A 2 \og(L/6t). We regard 
the variance Var(uist) as the quantity representing the intensity of the market-wide 
collective behavior, i.e. herding of market participants. 

In the MRW model, ugt is a stationary process. Therefore, the parameters A 
and L are constant against time. However, the intensity of the market-wide collec- 
tive behavior is thought to undergo a change depending on the market phase. We 
therefore evaluate the temporal behavior of the variance Var(cost) in a sliding time 
window [t — AT, t] with given width of To obtain the variance Var(uist) in each 

window, we must estimate parameters A and L. From analytical results of MRW 
model, the covariance of the logarithm of absolute return follows the equationp 1 



Cov(log(5X S t[i]),log(5Xst[i + k])) 



-A 2 log(|fc|/L) for 5t<.\k\<L 
L<\k\. 



(3-1) 



In Figure [3j an example of the covariance function for a time window is shown. The 
covariance function is fitted by the prediction of the MRW model (3T). Parameters 



A and L are estimated using this equation. We set width AT a fifth of the whole 
period, namely, AT = 39698 min, which is in the range L < AT < T. 
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Fig. 3. Example of covariance function for a time window. The semi-logarithmic plot of 
Cov(log(5Xst[0]), \og( 5Xgt [fc])) against lag k is shown. The straight line is the prediction of 
the MRW model p46. The estimated values are A 2 = 0.018 and L = 12975.43. 



In Figure |4j we show the temporal evolution of the variance A 2 log(L/5t). In 
general, for a non-stationary process, de-trending of the time-series must be pre- 
scribed. However, the de-trending method of is not determined uniquely. Therefore, 



6 



J. Maskawa 



we compare the result without de-trending with that using a thorough de-trending 
prescription^] In both cases, the sharp rises of the variance X 2 \og(L/5t), which 
mean the upwelling of the collective behavior of stock prices, are observed before the 
market crashes of Jan. and Oct. 2008. The periods of high volatility lasted several 
months after the crashes. This phenomenon is not observed for steep declines of 
Mar. and Sep. 2008, which were caused by the buyouts of Bear Stearns by J. P. 
Morgan on 17 March and the Lehman shock on 15 September, respectively. 
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Fig. 4. Temporal evolution of the estimated variance A 2 log(L/5t). Two results are shown. One 
is the result with no de-trending prescription and St = 1 min (black solid line). The other is 
the result with the thorough prescription described in footnote and St = 8 min (dashed line). 
The daily frequency of large 8M(t) values with amplitudes larger than the twofold standard 
deviation during the period (black lines) is also shown (gray line). 



As the final empirical result, we demonstrate that the sharp rise of the variance 
before the crash of Oct. 2008 is associated with the news which include the words 
"financial crisis." They have no severe impact on stock return alone but contribute 
to the pessimistic mood among stock market participants. The amount of such 
news increased after the Lehman shock preceding the market crash (Figure [5]) . The 
dynamics of the amount of news with specific words before and after the crash is well 
described by the precursory and the relaxation dynamics of a social system, which 
has been studied i r0 using time-series of daily views of video on YouTube. News 



Here we remove the local trend every 8 min as the thorough de-trending prescription. 
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including the words "financial crisis" increased toward the peak, taking about three 
weeks and showing slow relaxation lasting over the end of 2009, which is characteristic 
of an "endogenous" burst. In contrast, the news with the words "Lehman Brothers" 
show a sudden peak and rapid relaxation as an typical profile of " exogenous" burst E3 
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Fig. 5. Temporal evolution of the daily amounts of news which include the words "financial crisis" 
of the period from 1 Sep. to 31 Dec. 2008 (black solid line), the news with the words "Lehman 
Brothers" of the period from 1 Jul. to 31 Dec. 2008 (dashed line)0 The news sources are 
confined to Reuters news. The amounts of news strongly depends on the day of the week. To 
remove that effect, we divided the daily amount by the average value of the corresponding day 
of the week and multiply it by the average during the period. The temporal evolution of the 
variance A 2 log(L/St) with no de-trending prescription is also shown (gray line). Inset: Actual 
monthly count of both news of the period from Sep. 2008 to Dec. 2009. 

Figure [6] shows the increase of the variance Var(co$t) = X 2 log(L/5t) against the 
cumulative amount of news N n for the period from 16 Sep. to 8 Oct., which is 
expressed by the power-law 



Var{io 5t )(t) ~N n {t) a ", (3-2) 
where the estimated exponent a n = 0.19 ± 0.04 in this case. 
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Fig. 6. Growth of the variance Var(uist) = A 2 log(L/<5t) against the cumulative amount of news 
iVnfor the period from 16 Sep. to 8 Oct (circle). The solid line represents the best fit by eq. 
l ]3-2[ ). The estimated exponent a n is 0.19 ± 0.04 



§4. Conclusions 

As described in this paper, we have shown through an empirical study of the 
stock returns of the constituent issues of FTSE 100 index listed on London Stock 
Exchange for the period from May 2007 through Jan. 2009 that precursors of the 
market crash of Jan. and Oct. 2008 exist. A sharp rise preceding the crashes in the 
collective behavior of stock prices was measured by the variance of the logarithm 
of return, which is a parameter of the MRW model. This phenomenon signifies the 
upwelling of the market-wide collective behavior before the crash, which might reflect 
a herding of market participants. The crash of Oct. 2008 was synchronized with the 
increase of news including the words "financial crisis." The variance increases along 
with the cumulative amount of news according to a power law. 

It has remained unknown how to predict the time and the magnitude of steep 
decline such as that of a market crash. However, it might be accomplished through 
the accumulation of understanding about complex systems. 
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Appendix A 

Multifractal Random Walk (MRW) model 

Multifractal Random Walk (MRW) X{t) is a continuous process defined by the 
limit of the discrete random process X§t- 



X{t)= lim Xst(KstSt). (A-l) 

St— >0,t=K st 6t 

The discrete process Xgt(Kgtdt) is a stochastic volatility process that can be decom- 
posed into subprocesses 5Xg t : 

K st 

X St (K St 6t) = J2 Sx st[i\- (A-2) 

i=l 

The subprocess 6Xst(i) is described as 

SX 5t [i] = e st [i] exp(u S t[i]), ( A -3) 
where est is a stationary Gaussian white noise with variance o~ 2 5t and exp^^fz]) is 
the stochastic volatility. 

Bacry et al. show that the stochastic volatility process Xst(i) is a multifractal 
process p as 

M(q, St) = E(\5X 5t \ q ) ~ 5t(« (A-4) 
if ust[i] is a stationary Gaussian process such that E(oj$t) = —Var(ust) and 

Cov(ust[i],U8t\j]) = A 2 logj0 5t [K — i|], (A-5) 

where 

for \k\ < L/dt-1 



P5t[k] = )(\k\+m - (A . 6) 

II otherwise. 
Spectrum Q q is given by the formula 

C q = {q-q(q-2)X 2 )/2. (A- 7) 

To apply a MRW model to a time-series, three parameters o, A, and L must be 
estimated. 

Here, we check the applicability of MRW model to the market mode M(t) using 
the formulae of (A-4) and (A-7). Figure [7] shows the double logarithmic plot of the 
expectation values of q moments against the time scale 5t, and the spectrum £ ? , 
which is estimated by eq. (A-4). Those results mean that the MRW model is well 
applicable to the time series of the de-seasonalized market mode. 
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. 7. Multifractal properties of the time-series of the market mode, (a) Double logarithmic plot 
of the expectation values of q moments against the time scale St is shown. From top to bottom, 
the moment order q — 1,2,3,4,5. Time scale St varies from 1 min. to 4096 min. (b) Spectrum 
£ q (circle) and the best fit by the theoretical prediction (A- 7 1 (solid line) are shown. 
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